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The elementary excitations of a model Hamiltonian that 
captures the low energy behavior of a simple two-fold de- 
generate Hubbard Hamiltonian with Hund's rule coupling, is 
studied. The phase diagram in the mean-field limit and in 
a two-site approach reveals a rich variety of phases where 
both the orbital and the spin degrees of freedom are ordered. 
We show that besides usual spin waves (magnons) there ex- 
ist also orbital waves (orbitons) and, most interestingly, in 
a completely ferromagnetically coupled system, a combined 
spin-orbital excitation which can be visualized as a bound 
state of magnons and orbitons. For a completely degener- 
ate system the bound states are found to be the lowest lying 
elementary excitations, both in one- and two-dimensions. Fi- 
nally we extend our treatment to almost-degenerate systems. 
This can serve as an example that the elementary excitations 
in orbitally degenerate strongly correlated electron systems in 
general carry both spin and orbital character. 

I. INTRODUCTION 

The d and / wavefunctions of free atoms or ions are, 
besides Kramers (spin) degenerate, also 5 and 7-fold or- 
bital degenerate. In crystals this degeneracy may be 
lifted by the crystal field (interactions with the ligands). 
There are, however, interesting situations where this de- 
generacy is only partially lifted, so that in a high symme- 
try situation d- or /-levels remain degenerate. The best 
known examples are Cu 2+ (d 9 ), Mn 3+ (d 4 ) and Cr 2+ (d 4 ) 
in cubic symmetry (in octahedral surrounding) . Accord- 
ing to the famous Jahn- Teller theorenu this degeneracy 
should be lifted in the ground state. In a concentrated 
system this leads to the phenomenon known as,-jtiie co- 
operative Jahn- Teller effect, or orbital orderingHa. An 
interesting aspect of this phenomenon is the strong in- 
terplay between the orbital and spin (magnetic) order- 
ing. The orbital occupation determines the character 
of the magnetic exchange interaction (the Goodenough- 
Kanamori-Anderson rules)u; and, vice versa, modifica- 
tion of the magnetic structure, e.g. by an external field, 
may change orbital occupation and lead to a change in 
the crystal structuretJtj. 

The existence of orbital degrees of freedom, strongly 
interacting with the spins, does not only determine the 
orbital and spin structure in the ground state. It should 



also have important consequences for the elementary ex- 
citations of such systems. Thus, in addition to the collec- 
tive excitations of the magnetic subsystem -magnons, or 
spin waves-, also orbital excitations -orbitons= may ex- 
ist in this case. This was pointed out in refE and was 
studied for a specific model for the manganites in ref.B. 
Because of the intimate connection of spin and orbital 
degrees of freedom, one can also expect strong interac- 
tion and possible mixing of these two types of excitations. 
This problem was not addressed until now, and it is one 
of the aims of the present study. 

The systems with orbital degeneracy form a rather spe- 
cial class of compounds with very interesting and rich 
properties. Among them are, in particular, many com- 
pounds containing Cu 2+ (like actually the high T c super- 
conductors), manganites of the type Lai^Sr^MnOa with 
a colossal magneto-resistance, but also many other tran- 
sition metal compounds containing vanadium (LiV02, 
CaViOgpa, trivalent nickel (PrNi0 3 )B etc. The study 
of the elementary excitations in these compounds, be- 
sides of being interesting on itself, may shed some light 
on their unusual properties. 

In this work we consider the typical situation of mate- 
rials containing localized electrons with spin S= ^ with 
doubly degenerate orbitals. We can describe this situa- 
tion by the doubly degenerate Hubbard model: 

H = Ht + Hu + Hj (1) 

with 

Bt= £ *«4^ (2) 

<ij>,a,a,f3 

Hu= £ V afi n iaa n ll3a ,(\ - 8 a ,p8„ t „i) (3) 

Hj = £ -J^Sia • S (l - S a , fi ), (4) 

in which, besides the usual terms (electron hopping and 
on-site Coulomb repulsion), we also added the on-site 
(Hund) exchange interaction. 

In realistic situations the hopping matrix elements i™ 
depend on the type of orbital involved! - !. This leads 
to enormous technical complications. As here we are in- 
terested in the basic typical features of the excitation 
spectrum of our system, we treat a simplified symmetric 
model, assuming 
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*if = UjS a fi, (5) 

keeping only the nearest neighbor hopping for equal or- 
bitals. At the end of the paper we comment which mod- 
ifications of our results should occur for a more general 
choice of hopping integrals. In the case of strong in- 
teraction t <C U only the degrees of freedom of the lo- 
calized electrons with spin i and the orbitals are left. 
Consequently we can reduce the electronic model (Q) to 
a model describing coupled spins and orbitals. For the 
doubly degenerate case we can describe the orbital de- 
grees of freedom by an effective pseudo-spin T= |, so 
that e.g. an occupied orbital a corresponds to T z = \ 
and orbital j3 to T z = — i . The effective Hamiltonian has 
the generic form 

H = -J S J2 Si ' S, - J t £ T, • T, 

<ij> <ij> 

- U st J2 (S* • Si) (T, ■ Tj), 



(6) 



where for the model with t"^ given by (|^) the exchange 
constants have definite valuesB. Thus for this particular 
choice of hopping matrix elements the model (|l|) reduces 
to a double-spin Heisenberg model. 

The characteristic feature of our situation is the ex- 
istence and strong interplay of 'spins' S and T which 
emerges due to the third term of Hamiltonian (Q). As 
we shall see, this interaction leads to a significant mixing 
of spin and orbital interaction, giving rise, even, to the 
possibility of the formation of spin-orbiton bound states. 
Note that the form of the interaction between the spins 
and orbitals in (|l|) is a consequence of the Coulomb and 
Hunds rule interaction between electrons and is different 
from what one expects from simply coupling two Heisen- 
berg spin systems. As we want to study the generic fea- 
tures of the coupled spin-orbital system, we consider the 
general case, treating arbitrary values of the exchange 
parameters J s , Jt and J st . Notice that in the general 
situation with Hund's rule exchange included and with 
realistic values of the hopping integrals t"f the result- 
ing spin-orbital Hamiltonian of the .type (|3|) can contain 
terms anisotropic in the T-operatorscl and even terms lin- 
ear in T. However, even with the mare realistic hopping 
integrals there may exist situations^!] in which the sym- 
metry in orbital space remains similar to that of Hamil- 
tonian (||). For simplicity we consider such a symmetric 
double- Heisenberg model because, as it turns out, the 
conditions for the existence of the bound spin-orbital ex- 
citations are the most stringent exactly in this case, and 
if, as we will show, they exist in this case, it will be even 
more so with the T-anisotropy taken into account. 



II. GROUND STATE 

If our aim is to calculate the elementary excitations of 
a system described by the Hamiltonian (0) , we first have 



to know the ground state wavefunctions. We can easily 
calculate a phase diagram within two different approxi- 
mations: the mean-field and the two-site approximation. 
In the latter only two interacting spins and pseudo-spins 
are considered, so that the ground state is described in 
a way similar to a valence-bond state, which is a good 
starting point for low-dimensional spin systems. Start- 
ing from ground state wavefunctions obtained within one 
of these approximations, one can test the stability of the 
ground state by calculating its response to for instance 
spin waves and combined orbital- and spin-excitations. 



A. Mean Field Solution 

A common method used to gain some understanding of 
a Hamiltonian is a mean field approximation. In our case 
it is possible to separate spin and orbital degrees of free- 
dom by replacing operators we want to exclude from our 
consideration by the averages of their correlation func- 
tion that appears in the Hamiltonian. In this way we 
generate two mean-field Hamiltonians 



H MF = 



-Js ^ Si • Sj — 4J s t y < Tj • Tj > S,; • Sj 

<ij> <ij> 

-j's E s * • S; ( ? ) 



<y> 



and 



H t MF = 



~J t J2 T ■ T ~ X ' S <- H / ' T • T ./ 

<ij> <ij> 

- J t E T< •Tj- ( g ) 



<y> 




Jst<0 
Jt 1 



FIG. 1. Phase diagram of the model Hamiltonian in the 
mean field approximation. FF: both spin and pseudo-spin 
are ferromagnetically ordered; FA: spin ferro, pseudo-spin 
anti-ferro; AF: spin anti-ferro, pseudo-spin ferro; AA: both 
spin and pseudo-spin are anti-ferro ordered. J„ and Jt are in 
units of | Jst | - 

In Hamiltonian (Q) the orbital degrees of freedom are 
integrated out and in Hamiltonian (g) the average over 
the spin degrees of freedom is taken. In this way the spin 
and orbital degrees of freedom are decoupled. After doing 
the mean-field averagings we have, in some sense, thrown 
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away the interesting part of our problem and returned to 
a renormalized Heisenberg model, where the ground state 
properties and elementary excitations are well known. 
Nevertheless, as in a first approximation it is still a useful 
approach. The phase diagram in this approximation is 
given in figure |l|. 



B. Two-Site Solution 

In a mean-field approximation the short range interac- 
tions are averaged out, and if we want to go beyond this 
approximation and gain insight into local properties it is 
more useful to consider a few interacting particles. This 
is especially important here because the transformation 
properties of these terms of the Hamiltonian (|6|) are dif- 
ferent, and the more accurate account of the real singlet 
correlations (as compared to the mean field, essentially 
Ising-like, treatment) is essential. To this end case we 
consider two spins and two pseudo-spins and obtain the 
Hamiltonian 



-Hi, 2 — — Js&i ■ S2 — JtTi • T2 
- 4J rt (Si • S 2 ) (Ti-T 2 ) 



(9) 



This simple case can be treated exactly, and using these 
results we obtain the modified phase diagram. The com- 
binations of (pscudo-)spin operators can either be singlet 
or triplet and the ground state energy of the various com- 
binations correspond to the different phases in figure ^, 
for J st > and J s t < 0, respectively. The character of 
the first excited state for the various phases for J st > 
is shown in figure |^. A similar picture can be made for 
J s t < 0. In spite of the relative simplicity of the Hamil- 
tonian, there is a rich variety of ground- and first-excited 
states, even in the two-site approach. 
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FIG. 2. Phase diagram of the two-site solution of the model 
Hamiltonian. In the 11 phase both spin and pseudo-spin are in 
a triplet state, in the phase 10 the spin is triplet, pseudo-spin 
singlet, in the phase 01 the spin is singlet, pseudo-spin triplet 
and in the phase 00 both spin and pseudo-spin are in a singlet 
state. J s and Jt are in units of \J a t\- 

In the white, central part of this figure the ground state 
is singlet, both in the spins and orbitals. The first excited 
state, however, is not just the change of one singlet into a 
triplet, but corresponds to the state where both the spins 



and orbitals are in a triplet configuration. This is due to 
the interaction between the spin and orbital degrees of 
freedom and this very simple example shows that some 
of the excitations of a system with such an interaction, 
in this case the lowest one, carry both spin and orbital 
character. 

Jst>0 



□ 11 

□ 00 
Js Hoi 

□ 10 



FIG. 3. Character of the first excited state in the two-site 
system. J s and Jt are in units of \J s t\- 



III. FERRO-FERRO SYSTEM 

When J s , J t and J s t are chosen in a range where 
both the spins and pseudo-spins order ferromagnetically, 
a ferro-ferro phase, we can obtain exact analytical expres- 
sions for the elementary excitations. In a ferromagnetic 
Heisenberg model the exact ground state is the state with 
all spins pointing in the same direction, as opposed to a 
Heisenberg Hamiltonian with an antiferromagnetic ex- 
change where quantum fluctuations affect the Neel spin 
order. 

The coupling between the spins and pseudo-spins can 
lead to the formation of bound states of spin and pseudo- 
spin excitations. This is illustrated for the ferro-ferro sys- 
tem in figure [|. Let us consider Ising spins (and pseudo- 
spins) and examine the energy of an excitation of one 
spin plus one pseudo-spin. If the two excited spins are 
far away from each other, the part of the excitation en- 
ergy due to the spin pseudo-spin interaction (the third 
term in (||)) is 8| J s t|. When the excited spins are close to 
each other, on the nearest neighbor sites or on the same 
site, this excitation energy is 4| J st \ and 0, respectively. 
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FIG. 4. Excitation energy for Ising spins and pseudo-spins 
for the ferro-ferro system. When the spin and orbital are close 
together, the excitation energy is lower than in the case they 
are far apart. 
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This means that a magnon and an orbiton (orbital 
wave) have an attractive interaction and that in prin- 
ciple, depending on dimensionality and strength of the 
attraction, the magnon and orbiton can form a bound 
state. 



A. Equations of motion 

In order to establish whether it is indeed possible to 
obtain bound states that are combinations of orbital- and 
spin waves in the excitation spectrum of the model sys- 
tem, we use the equation of motion method, which pro- 
vides a simple, and in the ferromagnetic case exact, way 
to calculate bound state energies. Before addressing this 
issue, let us first examine a single (pseudo-)spin excita- 
tions. 

Starting from Hamiltonian (ra) and with a ground state 
|0 > with energy E where all (pseudo-) spins are aligned 
in the positive z direction, a single spin is excited (the 
derivation for single orbital excitations is equivalent). 
The equation of motion for this excitation is: 

(H - E )S m |0 >= u> s S- |0 >= [H, S-] |0 > . (10) 

The evaluation of the commutator and transformation to 
Fourier-space yields the following Goldstone modes: 

w,(Q) =2(J S + J s4 )^(l-cosQ-a) (11) 

a 

w*(Q) =2(J t + J sf )^(l-cosQ-a), (12) 

a 

where the lattice vectors are denoted by a. The presence 
of the coupling between the spin and orbital degree of 
freedom, parameterized by J s t , thus merely renormalizes 
the spin/orbital wave spectrum. 

The equation of motion for a combined spin and 
pseudo-spin excitation S^T~ is 

(H - E )S~T- |0 > = uS-T- |0 > 

= [H,S-T-}\0>. (13) 

We find for the Fourier transform of the equation of mo- 
tion 

{co- 7 (Q,q})A(Q,q) = 

-8J st ( — ) d V(cosQ • a/2 - cosq • a) 
2ir — ' 



J <ik(cos Q • a/2 - cos k • a)A(Q, k), 



(14) 



with 



A(Q,q) = S Q/2 _ q T Q/2+q |0> (15) 
7 (Q,q)=u; s (Q/2 + q)+u; t (Q/2-q) (16) 



where d is the dimensionality of the system. The total 
momentum of the excitation is Q and the relative mo- 
mentum of the combined spin and pseudo-spin excitation 
is k. In order to check if there is a self-consistent solution 
of the equation of motion, equation ( pT[ ) is iterated once. 
Then a set of equations is found that can be represented 
in a matrix M Qji g| of order d for hyper-cubic lattices. 
In a three-dimensional system a,/3 = x,y, z. The set of 
equations has a solution if, and only if, 



Det IS, 



M a>p \ = 



(17) 



with 



M a ,p{Q) = -8J st ( — ) d 
Zn 



(cosQ Q /2 - cos/c Q )(cosQ/3/2 - coskp) 
« k ,r\ i s 



If the system is one-dimensional the condition above re- 
duces to 



1 = 



8 J s t 
~~2n 



dk 



(cos Q/2 — cos fc) 2 



u-e.(k + Q/2)-£ t (k-Q/2) 



, (19) 



where the lattice spacing is set to unity. This expression 
resembles the condition for the existence of a two-magnan 
bound state in a pure ferromagnetic Heisenberg modeO. 



B. Bound states 

We restrict ourselves to the ferro-ferro phase, i.e. ac- 
cording to figure ||, 2J st < J s + Jt and J st > 0, and 
determine the bound state energies from condition ( |T7j ) 
and (|l^). The integrals appearing in the latter equation 
are, as in the two-magnon problem, by no means trivial 
and can only be determined analytically in special cases. 
In general these equations can be solved numerically. 
Before examining the integrals in detail, let us first turn 
to the excitations that lie in the continuum. The con- 
tinuum of excitations starts at the energy w c (Q) where 
the denominator of the matrix elements M a ^ starts to 
diverge. The condition for this is: 



w c (Q) = Min k [e s (Q/2 + k) + e t (Q/2 - k)] 



(20) 



The single (pseudo-)spin wave excitation energy always 
lies in the two-particle continuum: 

w a (Q) > w c (Q) 
w*(Q) > w c (Q). 

The explanation for this is that in a combined spin- and 
orbital-excitation the momentum is distributed over the 
two sub-systems and because the excitation spectrum in 
the ferro-ferro case is non-linear and concave ([1 — cos]- 
like), the energy of two excitations with smaller momenta 
can be lower than the energy of one excitation with large 
momentum. 
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First we consider the one- dimensional case. The inte- 
gral in equation ( |l9| ) can be reduced to the form 



(a — cosfe) 2 



1 + 6 cos k + c sin k 



dk, 



(21) 



which is known analytically. The condition that the in- 
tegral be equal to unity is given by the roots of a third 
order polynomial. Not so much is gained following this 
procedure, as the resulting expressions are very long and 
complicated. So let us consider some special momenta in 
the Brillouin zone. 

Q = 

For Q = the condition (|l^) for a bound state at u = 
reduces to: 2J st — J s + J t , which is exactly at the phase 
boundary of the ferro-ferro phase in the two-site phase 
diagram, figure ||. A bound state at negative u> appears 
when 2J st > J s + Jt- This simply means that the ground 
state is not stable (by exciting it energy is gained) , which 
can be expected as the two-site phase diagram shows that 
the antiferro-antiferro phase is the ground state in this 
parameter range. The ferro-ferro ground state is found 
to be exactly stable in the parameter ranges shown in 
figure ^, indicating that the two-site approximation gives 
a good prediction for the ground state spin and orbital 
order, whereas the mean-field solution fails to predict the 
right ordering. 



J.=J,=J„ 



S-O Boundstate 

S-O Continuum 
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0.671 
— >Q 

FIG. 5. Dispersion of the spin-orbital bound state, 
spin-orbital continuum, single spin and single orbital ex- 
citations in a one-dimensional ferro-ferro system where 
J s = J t = Jst- The single spin and single orbital excitation 
energies coincide in this case. The unit of energy is J = J s . 

Q = 7T 

By treating this special case, we can prove that there is 
always a bound state in the ferro-ferro system, at least in 
one-dimension. When Q a = ix for all a, at the corner of 
the Brillouin zone, the equations simplify considerably. 
The off-diagonal matrix elements in equation ( |l7| ) vanish 
so that Mq.,/3 = D6 a ,p. This yields 
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(tt) = 2{J S + J t ) - (J s - Jtf/^Jst (22) 
iol d {Ti) = 4( Min[J s , J t ] + J st ), (23) 



From these equations follows that (tt) < w)r (tt) for 
any J s , J t and J st , except when J st = \ J s — Jt\, where 
w b d ( 7r ) = ^c^ 71 ")- This proves that at Q — tt there is 
always a bound state. 

At uj(Q) -> uj c (Q) in ID 
For energies approaching the continuum (w(Q) — > 
w c(Q)), the integrand in equation (|l9|), diverges as fc -1 
(except for Q = 0), making the integral logarithmically 
divergent. We can conclude from this that for any Q 
(except for Q = 0) there is always a point between 
— oo < to < lo c {Q) where the integral is equal to unity. 

From the considerations above we can conclude that 
for the one-dimensional system in the range < Q < tt 
a spin-orbital bound state always exists and that this 
bound state is the lowest lying elementary excitation of 
the ferro-ferro phase of model (|6|) . It is by definition lower 
in energy than the spin-orbital continuum, which in turn 
is lower than the single (pseudo-)spin excitations. Before 
illustrating the statement above with numerical exam- 
ples, let us consider one more special case, namely: 

Js = Jt = Jst 

For these parameters the system is exactly at the phase 
boundary of the ferro-ferro ground state. Now equa- 
tion (^|) takes a particularly simple form, and the solu- 
tion for the bound state for the one-dimensional system 
is 



|(Q) = 2J st (l-cosQ). 



(24) 



For the lower bound of the continuum, the single-spin 
and the single-orbital excitations one obtains in this case 



uj* c (Q) = 8J st (l-cosQ/2), 
w s *(Q) = w t *(Q) = 4J st (l - cosQ). 



(25) 
(26) 



In figure]^ the dispersion relations of these excitations are 
shown. The spin-orbital bound state is always the lowest 
energy excitation of the system. Note that for small Q 
the spin-orbital continuum and the bound state are very 
close in energy and their energy difference is of the order 



of 4f Q 4 . 



J,=l= 2 l 



J,=0.5J=2J i: 




0.4ll 0.6JI 0.8n 71 0.271 0.471 0.671 
->Q ->Q 

FIG. 6. Dispersion of the spin-orbital bound state, 
spin-orbital continuum, single spin and single orbital exci- 
tations in a one-dimensional ferro-ferro system. The unit of 
energy is J = J s . 
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The numerical solution of the bound state equation 
for a one-dimensional system is shown in figure [| for two 
parameter sets. For small momenta the bound state is al- 
ways very close to the two-particle continuum. When J st 
is reduced, the single spin, single orbital and continuum 
shift down in energy, approaching the bound state energy. 
This can be expected since in the case when J s t = and 
J s = Jt the lower bound of the continuum and the single 
spin and orbital excitation spectra all coincide. 

In the right part of figure | a case is shown where the 
single spin and single orbital excitations have a different 
dispersion, i.e. J s ^ J t . The single orbital excitation 
is shifted up in energy with respect to the single spin 
excitation. The bound state, continuum and magnon 
excitations are close in energy, and in the limit that Jt — » 
oo all three coincide, as can be expected. 

In figure [?] a typical result for a two dimensional system 
is shown. It is found numerically that there always exists 
at least one bound state, also in two dimensions. Similar 
to the ID case, the bound state is only well separated 
from the continuum at the edge of the Brillouin-zone. If 
J s = J t it can be shown that 



uj b {n, 0) = 4(2 J s + 2 J st - ^j2J 2 st + 2J s J st + J s 2 ) (27) 
and 

Lu b (ir,ir) =4(2 J s + J st ), (28) 

so that at these points the bound state is considerably 
lower than the continuum. 




(0,0) 



FIG. 7. Dispersion of the spin-orbital bound state and 
spin-orbital continuum, single spin and single orbital exci- 
tations in a two dimensional ferro-ferro system. The unit of 
energy is J = J s . 



C. Almost Degeneracy 

In the sub-sections above we assumed that the two 
atomic levels are completely degenerate. Crystal fields 
generally split the two levels. Within the approach above 
it is fairly easy to treat this energy splitting. Suppose the 
energy difference between the orbitals a and is A. The 
Hamiltonian to be added to equation (fy) is 



H A = A/2^(n i!0%/3 - n i>at0l ) 

= A]T^, (29) 

i 

where the last equality follows from the definition of 
the T-operators. The level-splitting manifests itself in 
the pseudo-spin language as a magnetic field for the or- 
bitals. Carrying through the calculation for the excita- 
tions leads to renormalized spin-waves, orbital-waves and 
bound states: 



U! t = UJ t + A 
LO st = U>st + A. 

A/2 =/,=/=/„ 



(30) 
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FIG. 8. Dispersion of the spin-orbital bound state, 
spin-orbital continuum, single spin and single orbital exci- 
tations for an almost degenerate ID ferro-ferro system where 
A/2 = J s — J t — Jst- The unit of energy is J = J s . 

The spin- wave spectrum is not affected, but the "mag- 
netic field" for the orbitals causes a gap in the orbital 
excitation spectrum. This can be expected since due to 
the magnetic field there is no breaking of a continuous 
symmetry in the orbital case, and hence no Goldstone 
mode. This is reflected in the bound state energy be- 
ing gapped. The equalities above permit a convenient 
generalization of the results derived for the fully degen- 
erate system to the situation where the levels are non- 
(or almost-) degenerate, as illustrated in figure ||, where 
the dispersion of the bound state in the case of an orbital 
energy splitting of A = 2 J s is shown. 

From the results in this section one can also under- 
stand what will be the situation in the general case dis- 
cussed at the end of the introduction. Namely the orbital 
part of the effective spin pseudo-spin Hamiltonian is in 
the general case anisotropic, containing both terms of 
the type (^9|) and Ising-like terms TfT*. The situation 
then will be simular to the one discussed above: there 
will appear a gap in the orbiton spectrum (see alsofl). 
The spin-orbital bound state can still exist below the 
combined spin-orbital continuum (but in general above 
pure spin waves) and one can show that the conditions 
for their existence will be even less stringent than in the 
case of gapless orbitons. However, as can be seen from 
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figure ||, these bound states will not be the lowest exci- 
tations in the system, at least not in the whole Brillouin 
zone. Nevertheless they will lie lower then the orbital 
waves themselves, and they have definitely to be taken 
into account in a general treatment of properties of such 
systems. Note also that the orbital spectrum may still 
be gapless, even with realistic hopping integrals due 
to the orbital-lattice symmetry; this can lead to the ppj- 
ticular features of such systems like quantum disorderl!3; 
the nature of the elementary excitations plays a crucial 
role in this problerr£3. 



IV. DISCUSSION 

We studied the low energy excitations of a two-fold de- 
generate Hubbard model in the strong coupling limit with 
on average one electron per site and symmetric hopping 
integrals. We observe that besides the separate spin and 
orbital excitations, combined excitations (spin-orbiton 
bound states) can exist and can even be the lowest lying 
elementary excitation of the system. In general overlap 
integrals explicitely depend on the symmetry of the or- 
bitals; this, together with the Hund's rule coupling, may 
break the continuous rotational symmetry of the orbital 
channel. In this case one can expect that the eventual 
spin-orbital bound states, which are gappless in the sim- 
plest case, become gapped and may be the lowest lying 
excitation only at larger momenta. 

If not all interactions are ferromagnetic, there is a 
priori no reason to expect that the orbital and spin- 
orbital bound states behave qualitatively different from 
the ferro-case, although the quantum fluctuations might, 
or might not, destroy long range order£3ii3uJ. This in- 
teresting aspect of the system with antiferromagnetic in- 
teractions still deserves further study. 

The low energy collective modes of the orbitally de- 
generate Hubbard model certainly contribute to the ther- 
modynamic properties of the system, and should be ob- 
servable in for instance susceptibility and specific heat 
measurements. It should be stressed, however, that the 
elementary excitations with predominantly orbital char- 
acter are in general gapped, and therefore lead only to 
moderate changes in thermodynamic quantities. Exper- 
iments that are sensitive to higher energy scales might 
give direct evidence for the existence of orbitons and 
spin-orbital bound states. Orbital excitations, however, 
are strongly coupled to phonons and it might be difficult 
to distinguish between these two contributions in exper- 
iment. One of the possibilities to pin them down may 
be connected with the possible anomalies of the phonon 
dispersion relations induced by the orbital degrees of free- 
dom. Another possibility might be that, since the orbital 
excitations also couple to the spin degrees of freedom, 
variations in the spin order, for instance induced by an 
external magnetic field, will be reflected in the energy 
and dispersion of the orbital related excitations and con- 



sequently in the properties of the phonons mixed with 
orbitons. 

In conclusion, we studied a model Hamiltonian that 
captures the low energy behavior of a two-fold degen- 
erate Hubbard Hamiltonian. We presented the phase 
diagram in the mean-field limit and in a two-site ap- 
proach, revealing a rich variety of phases where both the 
orbital and spin degrees of freedom can be ordered (anti- 
)ferromagnetically. We have shown that in this case there 
may exist, besides usual spin waves (magnons) also or- 
bital waves (orbitons) and, most interestingly, the com- 
bined spin-orbital excitation which can be visualized as 
bound states of magnons and orbitons. In a fully degen- 
erate system the bound states are found to be the low- 
est lying elementary excitations, both in one- and two- 
dimensions. This shows that the elementary excitations 
in orbitally degenerate strongly correlated electron sys- 
tems in general may carry both spin and orbital charac- 
ter. 
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